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ABSTRACT 


Repeated or closely packed modal frequencies are common physical 
occurrences for vibrating structures which are complex or possess 
multi-planes of symmetry. The computation of the sensitivity to 
structural modifications for these frequencies and mode shapes is 
made difficult by the fact that the mode shapes are not unique, since 
any linear combination of eigenvectors corresponding to a repeated 
eigenvalue is also an eigenvector. 

This paper extends the work of Chen and Pan [1], who used modal 
expansion techniques for accommodating the sensitivity analysis of 
structures with repeated eigenvalues. Starting with a discussion of 
the physical significance of sensitivity analysis for repeated 
frequency modes, the paper presents a derivation of the governing 
equations for the derivatives of a repeated eigenvalue. This is 
followed with a small example to illustrate the results. An efficient 
computation procedure, based upon an expansion of Nelson's ideas [2] 
for large banded systems, is then proposed for systems with repeated 
or closely spaced eigenvalues. 
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IMPORTANCE OF THE PROBLEM 


The importance of obtaining gradients for eigenvalue problems 
stems from the fact that gradients, or derivatives with respect to 
system parameters, represent solution sensitivities. A knowledge of 
these sensitivities permits efficient design modifications, yields 
insight into the reasons for discrepancies between structural analyses 
and dynamic tests, and suggests model changes to improve correlations. 


Knowledge of Gradients: 

Yields Insight re. Parameter Sensitivities 


Permits Efficient Design Modifications 


Understand Test/Analysis Discrepancies 
Suggests Model Changes to Improve Correlation 
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WHEN DO REPEATED FREQUENCIES OCCUR? 


While a procedure for obtaining gradients efficiently was 
presented a decade ago by Nelson [2] , the problems associated with 
repeated roots have not been adequately addressed. 

The problem of repeated frequencies, or identical frequencies 
with different mode shapes, occurs in many physical situations. The 
most common circumstances under which multiple eigenvalues occur in 
engineering are cases where system symmetry exists, such as structures 
with two or more planes of reflective or cyclic symmetry (see Figure 
1) or axis symmetry (see Figure 2). 

It is also possible for repeated or closely spaced eigenvalues 
to occur when physical symmetry is not present, such as with classical 
wing flutter when the first bending and twisting frequencies coalesce. 


Coincidental Parameters (e.g, King Twist/Bending Flutter) 


Symmetry: Reflective, Cyclic, Axi symmetry 



Symmetrically Supported Mass 


Right Circular Cylindrical Shell 
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TECHNICAL BACKGROUND 


Assume [A] and [B] are symmetric n x n matrices and A* 
repeated eigenvalue with m+1 distinct orthogonal eigenvectors. 

{Z • ] is also an eigenvector corresponding to A ' vvhere 
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Symmetric Eigenvalue Problem 
Orthonormalization 
Multiple Eigenvalue 
Corresponding Eigenvectors 

Nonuniqueness of Eigenvectors 


( [A] - ) {x.} = {0} 

W T «1 {x.} - S t} 


A t * Repeats m + 1 Times 


l X t) ' ( x 4+1)-’ • ■ ■ • 'l x <:+ m| 


M 

- 2 °G( X J -w{ 


J--0 


M - 

" 1 1 

X - X ; 1 • ■ ■ 

1 ^ 

11 X [ + M 

lAJ - 

* c | 1 + 1 1 1 ' 

1 J 


is a 
Then 


l f 


200 


FRICTIONLESS PARTICLE IN A SHALLOW ELLIPTIC DISH 
PHYSICAL INTERPRETATION 


A simple physical interpretation of repeated eigenvalues was 
presented by Crandall [3 ] in which he considered a shallow elliptical 
bowl in which a frictionless mass particle is allowed to slide in the 
bottom of the bowl (left figure). The eigenvalue problem for this 
system consists in determining the paths and frequencies of 
back-and-f orth motion in which each motion is repeated on the same 
path. The natural mode solution is obviously along the major and 
minor axes of the ellipse. 

Next imagine that the elliptical bowl is gradually transformed 
into a spherical bowl (right figure). The eigenvalues will approach 
one another and any straight path, through the bottom of the bowl, is 
equally a natural mode. Thus, when m+1 eigenvalues coalesce, there is 
an infinity of mode shapes composed of a linear combination of the m+1 
dependent, but somewhat arbitrary, basis modes. 


Minor Axis 
Mode Path 



Elliptic Dish - Unique Modes 

Natural Frequencies Coalesce 
Distinct Mode Shapes Vanish 


Any Diametral Path Is A 
N atural Mode 



Spherical Dish - Non unique Modes 


As D ish Becomes Spherical 

As Ellipticity Disappears 
Crandall (ref. 3) 
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MODAL GRADIENT EQUATIONS PRESENT AN ENIGMA 


The modal sensitivity equations for a small change in a typical 
parameter, R, upon which certain matrix elements of [A] and [B] depend, 
are well known and summarized below. The problems are that they 
cannot be easily interpreted for the repeated eigenvalue problem since 
} is not unique, and matrix ([A] - [B] ) is not of order n-1 

but lower (i.e., n-(m+l)) depending on the multiplicity, m, of 

eigenvalue . 

( )'. 

D R 

x/ = {x ; J T ( [A] ' - MB] ' ) (x,) 

( [A] - A { [B] ) [x/J = {f-J 

{x ; j T [B] (x/j - b 

(fJ = - ( [A] [B] ) ' (x.| 



INTERPRETATION PROBLEMS 

There are ambiguities associated with the gradient equations 
since [x ^ f is not unique and A. depends upon which{X-} is chosen. In 
addition, the rank of ([A] - [B] ) is not n-1, but lower. 

Therefore, inclusion of the derivative of the normalization condition 
alone is not sufficient to uniquely determine^X.'^. 


Which (xj Should be Used? Is{ xJDifferentiable in R? 


Use of DifferentIx.) will Yield Different Results. 

L 


Rank of ([a] - \[b])is Too Low to Uniquely Determine |x/ 


Which Additional Equations Should be Used? 
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REPEATED FREQUENCY SENSITIVITY EQUATION 

To determine how the eigenvectors are perturbed by the 
infinitesimal change in R, we postulate an aribtrary vector [2 • } 
which is linearly composed of all the [xy } ( / = «* , * +1 , ... , 
i+m) and premultiply the eigenvector gradient equations by the 
transpose of all the eigenvectors corresponding to A; . 

This yields an auxiliary matrix eigenvalue equation in , 

which is of order m+1, the solution of which defines the specific 
eigenvectors, through the eigenvectors f )\, affected by the 

change in parameter R. 
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PHYSICAL INTERPRETATION 


Solution of the (m+1) eigenvalues and eigenvectors of the£D] 
matrix will yield the m+1 gradients of Ajas well as the eigenvectors 
{z to which they correspond. The figure below displays how the 
eigenvalues coalesce for a particular value of the parameter R and 
also shows how they correspond to different gradients. In general, 
there will be as many derivatives as there are curves intersecting at 
a particular parameter value R. 


M+1 

Solutions 


A 


[d] a: {cx} 



Determination of X; and Corresponding {<*} 
Uniquely Determines Mode for Gradient Sought 


i 
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PROPOSED SOLUTION PROCEDURE: OVERVIEW 


The solution procedure proposed is an extended version of 
Nelson's method for non-repeated roots. The method maintains the 
original matrix bandwidth while eliminating m+1 equation redundancies 
in the original eigenvalue system. 

The equations to be eliminated are determined by examining each 
eigenvector which corresponds to the eigenvalue whose gradients are 
desired, and establishing which elements are the maximum for each 
vector. These then correspond to which m+1 rows of ([A] - [B]) 
should be considered as redundant. If the maximum elements of any two 
eigenvectors correspond to the same row, then it is necessary to go to 
the next smaller element until a set of m+1 equations for removal is 
obtained . 

Rather than eliminate these rows and upset the system 
bandedness, we propose to extend Nelson's ideas by zeroing out the 
corresponding element of [F- } and then solve for JV; }. 

Based Upon Maximum Elements of{z : },{z L+l }' • ■ •'{% l+m} 


Zero-Out m+1 Rows and Columns of ( [A] - A' [BJ ) 


Zero-Out m+1 Elements from {F; } , j = L j l +1 , . . I +m 


Solve: ( [Al - A l [B] ) {v ; } = { F ; } 


Note: ( LA] - \[B] ) Has Same Bandwidth as ( [A] - A t .[B] ) 
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SOLUTION PROCEDURE: AUGMENTATION OF EQUATION 


The process described on the previous page yields a solution 
vector jVj } with m+1 zeroes. To this we append the m+1 eigenvectors 
[ Z o ] with appropriate constants C This combination is then 

substituted into the derivatives of che m+1 orthonormalization 
equations and the (m) additional optional equations to uniquely 

determine the m+1 constants C . . 

J K 


Introduce m+1 Additional Equations: 

Cl z j} T [ B ] «Q' =° 

L 1 J = l* t+l» C +i » • • • 

c+^\ 

Let {Z ; '} ={Vi} + I C Jjr { Z x } 

f-i 

And { Zj .} t [ B]{ Z ;} ={Z k } t [b] {z/} (Optional) 

Then C,j = + {Z K } T [ b] {Zj } -{Zj} t [b] K} 
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SIMPLE EXAMPLE: BASIC DESCRIPTION 


As a very simple application of this procedure consider a 
weightless straight bar of length L with end masses supported by 
linear springs. As the spring stiffnesses and masses approach one 
another, so do the two system frequencies. Thus, depending upon the 
method by which the normal modes are obtained, the mode shapes may 
vary. For the mode shapes presented below, either both masses vibrate 


simultaneously up and down together 


l x.} , 


or in opposition, fx~}. 


SIMPLE EXAMPLE: NORMAL MODES 


If we follow the procedure outlined earlier, and compute the 
system frequency gradients with respect to changes in m2, we obtain 
normal modes JZ.} and f Z ^ • These modes are associated with motions 
for which only moves, and motions for which only nij moves. 
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SIMPLE EXAMPLE: EIGENVECTOR GRADIENTS 

Following the computation procedure outlined earlier, the 
eigenvector gradients, [Z'-J and $Z' 2 ], for changes in m 2 are the null 

vector and the value of ^ [Z 2 ?/^ m 2 shown below. 


( z i j “ c 1 1 { z 1} + c i2 { z 2} =| 0 ) ' i - e. { Z x } Not affected by 

{z 2 'l = C^z,}* c 22 (z 2 } = / 0 \ 


M 2 CHANGE 


■ E. { 
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USE CAUTION WHEN WORKING WITH THE TOTAL DIFFERENTIAL 



The Total Differential May Not Exist Even Though the Partial Derivatives Do, 
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CONCLUSIONS 


The coordinate system and mode shapes initially selected for 
this example gave little physical insight regarding how the initial 
system would decompose due to a change in m~. Yet, this example 
yields a simple demonstration of the insight to De gained by following 
the proposed procedure. Thus, it is seen that the proposed 
mathematical procedure automatically yields m+1 distinct gradients for 
repeated frequencies and m+1 distinct modes, without requiring user 
dependence. 

The computational efficiencies suggested by Nelson- 3 have been 
expanded. These include: maintaining system bandwidth and 
consideration of only the m+1 repeated root frequencies. 


Eigenvalue Gradients For Repeated Frequencies 
Generally Yield Multiple Distinct Values 


Efficient Computation of Eigenvector Gradients 
For Repeated Frequencies is Possible 

i.e, Bandwidth May be Maintained 

Modal Expansion is not Necessary 

But, Must Introduce Modal Orthogonality Conditions in 
Addition to Normalization Condition 

Exercise Caution When Using a Total Differential Which is a Combination of 
Partial Derivatives for Repeated Frequencies, 
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